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FOREWORD 


This report describes the methodology and data used in developing fatigue and fracture 
analysis routines which were integrated into an existing computer program which per- 
forins weight and cost analysis of transport aircraft. This work was sponsored by the 
NASA Langley Research Center. The contributions of Mr. C. C. Poe, Jr., NASA/ 

LRC contract monitor, are gratefully acknowledged. He provided the stress intensity 
factor coefficients for stiffened panels used in the crack growth and residual strength 
analysis routines. 

The work was accomplished by the Structural Analysis and Mass 'Properties groups of 
the Convair Division of General Dynamics, San Diego, California, under the supervision 
of J. E. Ashton, Director of Structures and Design. The Program Manager was C. J. 
Tanner. The structural synthesis module, APAS, which includes the added fatigue and 
fracture routines was developed by G. S. Kruse. ""Phil Thorndyke aided in developing 
fatigue data and analysis procedures. Work on the weight and cost computer program 
is continuing under the leadership of B. H. Omari. Significant contributions were 
made by T. F. Reed and A. R. Stone in computer programming to produce the resultant 
Vehicle Design Evaluation Program,. VDEP-II. 
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SUMMARY 


This report presents the results of a research and development study performed under 
NASA Contract NASI- 12506. The objective of the study was to develop a "Computer 
Program to Assess the Impact of Fatigue and Fracture Criteria on Weight and Cost of 
Transport Aircraft. " Its intended use is as a preliminary design analysis tool to en- 
able the user to rapidly perform trade-off studies involving fatigue, fracture, static 
strength, weight, and cost. 

The approach consisted of developing analysis subprograms for fatigue life, crack 
growth life, and residual strength; and linking these to a structural synthesis module 
which in turn was integrated into a "Computer Program to Perform Cost and Weight 
Analysis of Transport Aircraft" developed under Contract NAS1-11343. The fatigue 
routine utilizes a flight profile to calculate damage and is named PRODAM. The 
second routine uses the same flight profile to calculate crack growth and is called 
PROGRO. The third routine RESIDS, determines residual strength of flawed structures. 

The structural synthesis module, APAS, was originally developed under a Convair 
IRAD study in 1972. Beam type structures (i. e. , wing boxes and fuselages) are 
synthesized in a multi-station approach to obtain optimum design for static strength 
for several structural configurations. As a user option PRODAM, PROGRO, and 
RESIDS may then be called to check the design for fatigue and fracture criteria. A 
redesign loop is employed to augment the section until the input criteria are met. 

Under this study the part definition module of the cost and weight analysis program 
(NAS1-11343) was expanded to be compatible with the upgraded structural synthesis 
capability. 

The resultant Vehicle Design and Evaluation Program is named VDEP-II. It is an 
accurate and useful tool for estimating purposes at the preliminary design stage of 
airframe development. A sample case along with an explanation of program applica- 
tions and input preparation is presented in the User's Manual. Table 1 is a summary 
of the program functional capability and Figure 1 is a program block diagram. 
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Table 1. Summary of the Program Functional Capability 
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Aircraft Balance 
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Area Ruled Fuselage Geometry 
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Structural Synthesis 
Parts Definition and Weight 


Manufacturing Cost 
Material Cost 
Engineering Cost 
Tooling Cost 

Total Vehicle Program Cost 
Return-on-Investment 
Airline Route Analysis 
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“synthesis' 
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DESIGN DATA 
PLAN FOR MS 


MATERIAL PROPERTIES 
CONSTRUCTION TYPES 
LOADS 


MANUFACTURING COSTS 
FIRST UNIT COSTS 
TOTAL PROGRAM COSTS 



Figure 1. Vehicle Design and Evaluation Program 
(VDEP-II) Block Diagram 
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SECTION 1 


INTRODUCTION 


With the steadily rising cost of aircraft production and operation, and with the large 
number of materials and structural design concepts applicable to flight vehicles, it 
becomes increasingly important to be able to assess the impact of aircraft design alter- 
natives in terms of weight and cost. Advances in technology have produced components 
of increased specific strength, and hence, decreased weight, but at the expense of re- 
quiring increasingly exotic materials and fabrication complexities. The result has been 
an inversion in the typical cost/weight relationship. One way cf combating rising oper- 
ating cost is to design new aircraft for longer service lives. This combination of op- 
posing constraints requiring high structural efficiency on one hand and long service 
life on the other has focused the attention of airframe designers on fatigue and fracture 
considerations. 

Fatigue analysis has been applied in varying degrees of rationality in the design of most 
transport aircraft developed in the last three decades. Over-simplified approaches 
(such as setting an arbitrary limit on allowable tensile stress) have frequently been 
used in preliminary design. On the other extreme, methods have been devised, and 
used to solve special problems, which account for non-linear damage accumulation 
under complex load spectra. In this program an intermediate approach has been 
developed which permits a realistic evaluation of fatigue criteria and at the same time 
provides the necessary computational efficiency. This approach employs a simplified 
flight profile, component S-N data, and linear damage summation. 

A major obstacle to applying fatigue analysis in preliminary design is the sparsity of 
S-N data for the many structural configurations which need to be evaluated. Fatigue 
life is affected by a wide range of parameters that include cyclic stress, mean stress, 
product form and orientation, temperature, environment, structural geometry, notch 
effects, metallurgical effects, and surface finish. An almost infinite variety of materials 
and fabrication methods would have to be tested to build an adequate S-N data bank. Since 
this is out of the question, the fatigue life of most aerospace vehicles is verified by tests 
and analysis performed late in the development of the structure. This after-the-fact 
analysis cannot produce optimum structures. Tradeoff studies performed in the pre- 
liminary design phase are compromised by lack of rational fatigue considerations. Two 
methods for generating component S-N curves from unnotched coupon data have been 
applied in this program. Results are presented in Appendix A and B. A third approach 
utilizing an S-N equation directly within the computer program is proposed for future 
development. See '’Recommendations'*, page 5-1. 

Fracture mechanics is still a developing technology even though problems with slow 
crack growth as well as catastrophic cracking have afflicted various types of structure 
for many years. World War II military aircraft structures were designed using 
methods that accounted for the stiffness and static strength of the undamaged structure 
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only. Low-strength, highly ductile alloys such as 2024-T3 were used. Fabrication 
consisted primarily of forming and riveting sheet metal. Conventional sheet- stringer 
construction gave sufficiently stiff covers for the subsonic flight loads, along with 
reasonable structural efficiencies. The combination of low-strength ductile material, 
riveted sheet-stringer construction, , and thin sheet metal sections made these 
structures inherently highly damage tolerant. 

Post-war military aircraft were designed for supersonic flight, heavier payloads, 
and increased structural efficiency to save weight. Thick, monolithic, integrally 
stiffened, machined components replaced sheet- stringer construction to provide 
smooth, nonbuckling covers. High-strength materials such as the 7000 series alumi- 
num alloys and heat-treated steels replaced the low-strength ductile materials previously 
used and resulted in increased design stress levels. The combination of construction 
type, higher design stresses, more severe usage and more fracture-prone material 
produced an ever-increasing number of catastrophic service and test failures.*- In 
most instances ,, the failures originated at points of high stress concentration, or at 
metallurgical defects. Most often, the fatigue cracks, before failure, were hidden 
or were not large enough to be detected by close visual inspection. As a result, 
most high-performance military aircraft used during the 1950's and I960' s have 
gone through costly structural test, structural modification, and life monitoring pro- 
grams to remain flightworthy. Primary reliance has been placed on a safe-life 
approach to maintain serviceable aircraft. Extensive element, component, and com- 
plete airframe test programs are used to prove the safe-life of the structure. 

Commercial aircraft designed before and during World War II used material and 
structural concepts similar to those military aircraft. Again, the combination of 
low-strength ductile material and thin sheet metal sections made the designs highly 
damage tolerant. Later, two ill-fated airplane programs caused commercial air- 
craft designers to adopt a fail-safe design approach. These were the Martin 404, 
with heavy monolithic spar cap brittle fractures; and the De Haviland Comet pressur- 
ized cabin failures. Since these two instances, commercial aircraft have been de- 
signed with a strength capability of 80% of limit flight loads and 100% of normal cabin 
pressure with single structural members failed in compliance with FAR Part 25. These 
obvious partial failures are usually chosen to be sufficiently large to be apparent during 
normal preflight inspections. 

The design task of providing damage tolerance was made easier by the reduced stress 
levels and increased toughness provided by the almost exclusive use of a low-strength, 
ductile aluminum alloy (2024-T3) in these tension critical applications. The proof of 
the fail-safe load capability of the flawed structures was accomplished largely by 
fail-safe tests. 

In 1958, the Department of Defense, motivated by several fracture problems in rocket 
propellant tanks , initiated a study of fracture phenomena within the American Society 
for Testing and Materials (ASTM). The ASTM, through its E-24 Committee, provided 
considerable impetus to the advancement of fracture mechanics analysis. The com- 
mittee's work has resulted in standardized test methods and test specimens for the 
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proper measurement of material toughness (K Ic ). Reference 1 presents these results 
and has been used extensively as the text on fracture mechanics by the aerospace com- 
munity. Reference 2 is a comprehensive summary of fracture data for commonly 
used aerospace materials. 

Irwin (Ref. 3 ) developed equations that related the crack phenomenon to terms that are 
common to the designer or stress analyst (or, in his words, "suitable for stress analysis"). 
Later, Paul Kuhn developed practical guidelines that more or less guaranteed fail-safe 
structures without rigorous use of linear elastic fracture mechanics (Ref. 4 ). While 
Irwin's work provided insight into the cause of previous brittle fracture service prob- 
lems, his analysis methods were applicable only to the simplest structural components. 
Similarly, Kuhn's guidelines were not always compatible with the emphasis on improved 
structural efficiency and minimum weight. 

Stress intensity factor solutions for numerous load/ geometry combinations are now 
available. Reference 5 summarizes many of these. Special computer programs have 
been developed to solve additional complex geometries. Methods and programs are 
also available to predict crack growth under sustained and cyclic loading. The problem 
of crack propagation thru stiffened panels is of particular interest to aircraft designers. 
Methods and data required to solve this type problem are presented in Reference 6 
and are the basis of the crack growth and residual strength analysis routines developed 
in this program. 

There is a danger that this proliferation of design criteria, and analysis methods re- 
lating to fatigue and fracture criteria may reach a point of diminishing return, and 
even negative return. Therefore it is appropriate to evaluate and balance the possible 
improvements in service life and operating costs versus weight penalties and manu- 
facturing costs. That was the primary objective of this program — to develop a computer 
program to assess the impact of fatigue and fracture criteria on the weight and cost of 
transport aircraft. 

The starting point of this effort was a computer program developed under NASA Con- 
tract NASI -11343, "A Computer Program to Perform Cost and We ight Analysis of 
Transport Aircraft", (Reference 7 ). This program, intended for use at the prelimin- 
ary design level, incorporates both batch mode and interactive graphics run capability. 

The basis of the weight and cost estimation method developed is a unique way of pre- 
dicting the physical design of each detail part of a vehicle structure at a time when only 
configuration concept drawings are available. Encompassing the areas of manufacturing 
and material cost, engineering cost, tooling cost, total vehicle program cost, and re- 
turn-on-investment, this program represents a significant extension and refinement of 
the methods originally formulated in Reference 8. 

Weight data are generated in four areas of the program. Overall vehicle system weights 
are derived on a statistical basis as part of the vehicle sizing process. Theoretical 
weights, actual weights, and the weight of the raw material to be purchased are derived 
as part of the structural synthesis and part definition processes based on the computed 
part geometry. 
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The manufacturing cost analysis, based at the individual detail part level, is made by 
considering the actual manufacturing operations required to produce that part. A list 
of shop operations is called out with each detail part, and a series of equations associated 
with each operation is used to compute the shop hours necessary to make the part. 

By applying the appropriate labor rates to the calculated hours, the direct and indirect 
manufacturing labor costs are found. Material costs are computed based on the amount 
of material required to manufacture each part. 

Engineering costs are computed based on the number of manhours necessary to perform 
the various tasks associated with the development and production of aircraft. The com- 
putation has as its basis equations originally developed by Levenson and Barro of the 
Rand Corporation (Reference 9 ) . Initial engineering hours are broken down and dis- 
tributed among the various engineering disciplines based on studies made of historical 
data. 

Tooling costs are computed as a function of the number of basic tool manufacturing 
hours, initial and sustaining aircraft production rates, and tooling labor rates. Basic 
tool manufacturing hours are derived as a function of the number of dissimilar parts 
to be produced, the average number of tools required per dissimilar part, and the 
average number of hours required to produce each tool. 

Total vehicle program costs are computed based on a cost model that was assembled 
primarily from the work of Kenyon (Reference 10 ) . Cost elements that are computed 
elsewhere in the program are brought across and substituted into the model. A learn- 
ing curve approach is utilized to derive costs of a given unit or lot as a function of the 
first unit cost. ^ 

A comprehensive measure of the total economic viability for a commercial transport 
operation is reflected in the return-on-investment analysis. Direct operating costs 
are computed using the 1967 Air Transport Association formula updated to 1972 cost 
levels. Indirect operating costs and return-on-investment are computed by applying 
aircraft acquisition and direct operating costs to a defined traffic structure. Output 
includes direct operating costs, indirect operating costs, revenue, load factors, profit, 
return-on-investment, and fleet size. 

An important feature of the program is its capability to make trade studies from several 
levels of consideration. For example, weight and cost data can be related directly to 
key system parameters at the vehicle mission level such as payload, speed, range, and 
landing field length requirements. At the vehicle configuration level, data can be re- 
lated directly to surface areas, span, sweep, taper, etc. , and fuselage length, slender- 
ness, etc. At the major component level comparisons can be made between different 
materials, modes of construction, detail part make-up, etc. Tradeoffs can be made 
to determine the overall vehicle weight and cost sensitivities at each of these levels, 
and in this manner the proposed aircraft design may be further and further refined 
down to high degree of detail. Thus, engineering functions can gain insight into the 
cost effectiveness of alternate aircraft systems, perform design trade studies, as well 
as perform studies to determine the impact of fatigue and fracture criteria. 
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SECTION 2 


STRUCTURAL SYNTHESIS 

The structural synthesis procedure used by VDEP-H is known as APAS. This pro- 
cedure was originally developed earlier under a company sponsored IRAD program 
and is reported in Reference 11. 

The pertinent information of that report is restated in this section for convenience. 
Also reported in this section is a complete description of the technical approach 
used to include fatigue and fracture criteria in the structural synthesis procedure. 

2. 1 COMPONENT GEOMETRY 

The geometry of each component, (fuselage, wing, horizontal and vertical stabilizer) 
is represented by the coordinates of a set of nodes at each of the various stations 
along the component. This nodal geometry describes the shape of the component 
which is used for the computation of section properties. 

2. 1. 1 FUSELAGE NODAL GEOMETRY . The fuselage is represented by 18 nodes 
at each of 10 stations. The nodes are located at 20 degree intervals around the 
fuselage. Nodal geometry for a typical transport fuselage is presented in Figure 
2-1. Nodes are numbered starting at the top centerline and proceeding clockwise 
looking aft. 




Figure 2-1 Fuselage Nodal Geometry 
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Figure 2-2 Aerodynamic Surface Nodal Geometry 


2.1.2 AERODYNAMIC SURFACE NODAL GEOMETRY . The wing, horizontal and 
vertical stabilizer are represented by 10 nodes at each of 10 stations. The nodal 
geometry describes the box structure for a surface with up to 5 spars. The nodes 
are numbered beginning at the upper sparcap of the front spar and proceeding clock- 
wise to the lower front sparcap. A typical surface nodal geometry is presented in 
Figure 2-2. 

2.2 STRUCTURAL ELEMENTS, RIBS AND FRAMES 

2.2.1 STRUCTURAL ELEMENTS. Structural elements include skin panel, spar 
webs and spar caps. Each element is described by a type number and by from 
one to eight dimension variables. The dimension variables are of two types, thick- 
ness variables and non-thickness variables such as stiffener spacing, stiffener 
height, corrugation angle, etc. In general, non-thickness variables may have either 
equality or inequality constraints imposed, whereas thickness variables may have 
equality or lower bound inequality constraints imposed. 
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2.2. 1. 1 Skip Panel Elements. The structural synthesis program includes twelve 
types of panel elements as presented in Figure 2-3. The stiffeners on panel types 
one through nine are assumed to be oriented parallel to the elastic axis of the struc- 
ture. The 0 degree ply of panel type 12 is also assumed to be parallel to the elastic 
axis. 



Figure 2-3. Skin Panel Elements 
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2.2. 1.2 "Spar Web" Elements. The structural synthesis program contains seven 
types of "spar web" elements. Four of these are truss type elements, two are 
stiffened webs and the remaining one is a corrugated web. These elements are 
presented in Figure 2-4. "Spar Web" elements are assumed to resist only sliear 
and crushing loads, the axial stiffness of these elements is assumed to be zero for 
the purpose of computing section properties. 
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2.2. 1.3 Spar Cap Elements . Four types of spar caps are currently available. They 
include integral tee and angle and riveted tee and angle as shown in Figure 2-5. 




TYPE 1, RIVETED ANGLE TYPE 2, RIVETED TEE 

Figure 2-5 




TYPE 3, INTEGRAL ANGLE i TYPE 4. INTEGRAL TEE 

Spar Cap Elements 


2.2.2 RIBS. The types of ribs available within the program are presented in 
Figure 2-6 . The ribs are comprised of caps and webs or truss elements. Rib caps 
are sized to react a moment at the rear , spar due to the loading on the surface aft of the 
rear spar. Rib webs are sized to carry shear and to support crushing loads. 



built-up 

WEB 
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Figure 2-7 Typical Ring Frame 


2. 2. 3 FRAMES . A typical ring frame is shown in Figure 2-7. The frames 
are sized so that the outer flange clears all of the skin stiffeners. The inner flange 
is maintained at 14 cm (5.5 inches) from the outer skin contour. The frame is sized 
using Shanley's criteria to set a minimum frame bending stiffness. The frame is 
set to minimum gage for non-critical areas. 

2 

C f MD 

El = — — Shanley's criteria (Ref. 18) 

!j 

where: 

El = frame bending stiffness 

M = maximum resultant fuselage bending moment, yM 2 + M 2 

X z 

D = fuselage diameter 

L = frame spacing 

= fit coefficient (. 00025) 
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2.3 FLIGHT PROFILE AND LOAD SPECTRUM 


The fatigue load spectrum defines the number of times that incremental loads of given 
magnitudes are encountered during the design life of the aircraft. Experimental data 
is available which defines the probable magnitudes and frequency of occurrence of 
these incremental loads as a function of aircraft type, configuration parameters, and 
flight parameters. 

The configuration and flight parameters are defined using a typical flight profile, which 
is divided into segments. Parameters are averaged for each segment, and these aver- 
age values are used in finding the incremental loads. See Figures 2-8 and 2-9. 

The typical flight profile used for fatigue and flaw growth analysis is based on medium 
range operation of a contemporary transport aircraft. 

The parameter values for each segment are listed in Table 2-1 . The segments are 
divided into subsegments, with each subsegment representing a particular magnitude 
of incremental load. Using the segment parameters and the subsegment load, fre- 
quency of occurrence of the incremental load is found for each subsegment using the 
methods and information in Reference 13. 

For gust loads, curves showing gust velocity vs frequency of occurrence are found in 
Reference 13, Figure C13-32 through C13-37. From Reference 13, Page C13-24, 

Ag= mSV U K P o/2W 
e de g 

For speeds below critical Mach number: 

. 88 p 


g 5. 3 + p 

= incremental load factor 
= slope of lift curve 
= wing area 
= equivalent airspeed 

= derived gust velocity 

\ 

= gust alleviation factor 

= air density at sea level 

=* aircraft weight 

= aircraft mass ratio 

= acceleration of gravity 
c = mean geometric wing chord 

p = air density 


where 
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OCCURRENCES PER FLIGHTHOUR 


Figure 2-8. Typical Flight Profile 
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Figure 2-9. Typical Segment Load Frequency Curve 
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solving for 




Ag 


mS V K 

e g 


P o/2W 


U , is then calculated for each subsegment, and the curves of Figure C13-37 are used 
de 

to find the frequency of occurrence. 

For maneuver loads, Figure C13-41 of Reference 13 shows incremental load factor vs 
frequency of occurrence. For taxi loads. Figure C13-46 of Reference 13 shows incre- 
mental load factor vs frequency of occurrence. Incremental load factor vs frequency 
of occurrence for landing loads was averaged from data for two commercial transport 
aircraft. 

The resulting fatigue load spectrum is shown in Table 2-2. The number of cycles is 
based on 10,000 flights. The variation between cycles and flights is linear, so that 
linear ratioing of cycles and design life is valid. 


I 

Table 2-1. TYPICAL TRANSPORT FLIGHT PROFILE 


Segment Description 

Gross 

Weight 

(kiDsl 

Altitude 

CftxlO 3 ) 

* 

KEAS 

H 

§§§ 

1 Taxi; takeoff run, ldg. roll 

358.3 





2 Climb (Flaps down 25° ) 


0-5 


1 

8.15 

3 


■ 

■ 

5-10 

250 

.435 

12.51 

4 





10 - 20 

340 

.684 

47.53 

5 Climb 

352.2 

20 - 35 

319 

' 

150.70 

6 Cruise 

341.4 

35 

273 

.85 

395.58 

7 Descent 

335.4 

35 - 20 

319 

.836 

53.29 

8 



L_ 

20 - 10 

340 

^^9 


9 

1 ( Flaps down 15 ° ) 


M 

10-5 



SHIS 

10 Descent (Flaps down 50° ) 

335.4 

5-0 

IS 

.355 

ira 

11 Climb (Flans down 25° ) 


m 


igSl 

KSH 

12 



■ 

5-10 

250 

9 

i2. 51 

13: 

' 





10 - 20 

340 

.684 

47.53 

14 Climb 

352.2 

20 - 35 

319 

.836 

150.70 

15 Cruise 

341 

.4 

35 

273 

.85 

395.58 

16 Descent 

335.4 

35-20 

319 


53.29 

SSH 

- 



20-10 

340 

.684 

33.33 


(Flaps down 15° ) 



10-5 

250 

.435 

23.45 

19 Descent (Flaps down 50°) 

335.4 

5-0 

223 

.355 

14.11 

20 Landing 

334.7 

09 

'128 




♦Knots equivalent airspeed 
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Table 2-2. TYPICAL TRANSPORT FATIGUE SPECTRUM - CYCLES PER 10, 000 FLIGHTS 
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2 . 4 EXTERNAL LOADS 


Net limit loads due to the air loads, inertia loads and landing gear loads of various 
flight and ground conditions are input to the program. 



Figure 2-10. External Loads Sign Convention 

The loading conditions are separated into two groups. The first group consists of 
from one to six conditions. These conditions are specified by the user and are used to 
size the structure so as to preclude static strength failures and to meet residual strength 
requirements. The second group consists of the six conditions listed in Table 2-3. 

These conditions are used to define the fatigue stress spectrum described in Section 
2.7. 

Table 2-3. Fatigue Spectrum Loading Conditions 


Condition 

Number 

Description 

1 

1G Taxi 

2 

1G Level Flight 

3 

2G Vertical Gust 

4 

2G Maneuver 

5 

1G Landing Impact 

6 

Maximum Pressure 


(Fuselage) 


Each loading condition defines the six components of load (AX, XS, ZS, TOR, XM, 
ZM) at up to twenty stations along the structure. The sign convention used is pre- 
sented in Figure 2-10. A typical fuselage loading condition is illustrated in Figure 
2-11 . Steps in the loading curves are represented by repeating stations with the two 
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2 . 5 STRUCTURAL DESIGN PROCEDURE 


The structural design procedure starts with the input design, then through a series 
of design analysis and redesign iterations produces a final design. The iteration 
process continues until the lightest weight design which satisfies design criteria and 
constraints is found. 

The design analysis involves the comparison of applied stresses and allowable stresses. 
The internal load solution described ini Section 2. 7 is used to calculate the applied 
stresses. The box beam internal load solution was selected instead of a finite element 
solution in order to keep computer execution time at an acceptable level. This selec- 
tion restricts this procedure to relatively clean beam-like structures. However, de- 
coupling of internal loads from one station to the next is basic to the box beam theory. 
Hence, the overall design problem is reduced to a series of cross-section design 
problems at any number of desired locations along the structure. 

The procedure used to design the cross section is a two part procedure. The first 
part, the section sizing procedure, adds or subtracts material from the structural 
elements of the cross section in order to produce a zero margin or minimum gage 
design. The second part employs a non-linear programing technique to minimize the 
"criticality" of each element while maintaining a constant weight design. This element 
optimization procedure, is then iterated with the section sizing procedure until the 
design converges. Convergence occurs when two successive iterations produce a 
change in weight that is within a specified tolerance. 

2.5.1 SECTION SIZING PROCEDURE . During this part of the design process only 
the thickness variables are changed. Figure 2-12 illustrates a typical cross-section 
design problem. 
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The equivalent thickness of the panel, t , is defined as the panel section area divided 
by its width. For the panel illustrated, t is given by 

V~ *1,1 + b i,2 * ^,2^,1 = (Area of panel i)/ds. 

The following steps are used to find critical t 's. A critical t is defined as the t for 
which the panel's critical margin of safety is zero or for which the most critical 
margin of safety is positive and the panel is minimum gage. 

1. Determine the section properties of the cross-section. 

2. Determine the internal loads of each element of the cross section for each 
load condition. 

3. Calculate the applied stresses. 

4. Perform a stress analysis of each element for each load condition and produce 
the structural margins of safety. 

5. Starting with the 1st element of the cross-section and proceeding to the last 
find the critical constraint parameter, i. e. , lowest margin of safety or minimum 
gage margin. 

6. Change the t of each element to drive it towards its critical t. 

7. Return to step one and iterate the process until all panels have critical t's. 

The technique employed in step 6 to predict new t's is described below. The new t 
of a structural element is predicted bv passing a second degree parabola through the 
points (0, -1), (t k_1 , CPk - l) and (t*S ci^) then solving for the proper root of the 
resulting equation. The critical constraint parameter ( CP, is the margin of safety or 
the side constraint margin defined in Section 2.5.2. The superscript is used to denote 
the iteration number; i. e. , t^ is the current t , CP* 1 is the critical constraint parame- 
ter for t = t* 5 , and (t^ CP^" 1 ) are corresponding values for the previous iteration. 

The process, illustrated in Figure 2-13, is started by passing a straight line through 
the points (0, -1) and (t 1 , CP 1 ). 

2.5.2 ELEMENT OPTIMIZATION PROCEDURE . The object of the element optimiza- 
tion procedure is to adjust detail dimensions of an element so as to make the most 
efficient use of the material while maintaining a given weight. As an example refer 
again to Figure 2-12. The panel element shown contains four design variables, 

V V b l andb 2- 
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Figure 2-13. Section Sizing Procedure 


The object is to find the optimum set, i.e. , that set which represents the "least 
critical" panel possible. Since this procedure requires that t remain constant, there 
are only three independent design variables. Given any three variables the fourth 
is found by solving the following equation for the appropriate variable. 


l l + *2 b 2 /b l 


The "least critical" design is defined as that design for which the following criticality 
parameter is a minimum 


L J M 



where 

denotes the failure mode 
j denotes the loading condition 
m denotes the sub element of each panel element 

MS = margin of safety 

MC = side constraint margin e. g. , (t -t )/t . 

m min m mm m 
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x > C 


F(x) = ~ 
x 

i r x x t 

F(x) = r Lc ( e ■ 3 -) +3 -l 


where 


€ = 



A typical plot of the criticality function, F , is illustrated in Figure 2-14. 



Figure 2-14, Criticality Function 


There are a number of techniques available to minimize P, the method used by 
APAS is the Fletcher-Powell-Davidon; unconstrained minimization technique (Ref. 14). 
P is the objective function and [t , b^, b^] is the design variable vector for the 
example problem. Optimization is performed on each element of the cross section in 
order to minimize weight. 

2.6 STRUCTURAL ELEMENT SYMMETRY GROUPS 

A symmetry group is a group of structural elements which have identical designs. 

When a number of structural elements are placed into the same symmetry group only 
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one design is produced. The design of the element respects all of the margins of 
safety for all of the elements in the group. This technique provides a means by which 
fuselage centerplane symmetry can be respected without duplicating reversible loading 
conditions. It is often desirable to make adjacent panel elements identical for ease in 
manufacturing. This can be accomplished with symmetry grouping also. Since the 
use of symmetry groups reduces the number of independent design variables of the 
structure it can be of significant use in reducing execution time and should be employed 
wherever possible. 


2.7 STRUCTURAL ANALYSIS 

This section presents the techniques used to calculate the applied stresses, including 
the fatigue stress spectrum and the methods used to calculate margins of safety for 
static strength, fatigue, flaw growth, and residual strength criteria. 

2.7. 1 INTERNAL LOADS ANALYSIS. The internal loads analysis is based on 
classical box beam theory (Ref. 13). The assumptions made are; plane sections 
remain plane under the action of bending moments and axial loads, cross sections are 
free to warp when torque is applied, and the structure obeys a linear elastic stress- 
strain law. 


2. 7. 1. 1 Axial Stress. The axial stresses are made up of stresses due to axial loads 
and stresses due to bending moments. The equation used to calculate the axial 
stresses is, 

MI-MI MI-MI p 

a = -^LiE_^x (x _- ) + _z_xz x_E*( Z -E ) + f 

I I - 1 2 I I - I 2 

A XX ZZ XZ XX zz xz 


where, 


M = Net bending moment about a horizontal axis passing through 
the centroid 

M = Net bending moment about a vertical axis passing through the 
centroid 

P = Axial load 

x,z = Coordinates of the element, see Figure 2-15 

x, z = Coordinates of the centroid, see Figure 2-15 


I ,1 ,1 , A 

xz xx zz 


Section properties, see Section 2.7. 2. 
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2 . 7. 1.2 Shear Stresses. The shear stresses resist the external shear forces and the 
torque applied to the section. Under the basic assumptions the shear flow is calculated 
using a VQ/I distribution for the shear forces. The resultant applied torsion is due to 
the applied torque, TOR, and the couples resulting from shifting the shear forces 

XS and ZS to the shear center. This net torsion is resisted internally by a shear 
flow which is distributed according to a T/2A distribution. In the case of multiple 
cell structures, such as multi-spar wings, the cells are assumed to have equal twisting 
angles. For a further description of the method see Paragraph 17.9 thru 17.11 of 
Reference 19. 

2. 7.2 SECTION PROPERTIES. Section properties of the cross section of the wing 
or fuselage are calculated at each station where structural sizing is performed. 

These properties are used to calculate the internal loads distribution and to provide 
stiffness information. In order to simplify the calculation of section properties 

the following assumptions are made; ( 1 ) the material which resists bending moments 
is assumed to be smeared uniformly between nodal points, ( 2 ) only the skin and shear 
webs are effective for resisting shear loads and torsion. The following equations are 
used to calculate section properties , (see Figure 2-15). 

A 

x 

z 

I 

XX 

I 

zz 

I 

xz 
J 

where, 

x, z are the coordinates of the Incremental area da 
A is the total area of the cross-sectional material 
x,z are the coordinates of the centroid of A 


/da 

A / xda 
3 /zda 

f A " I 2 

J z da - A • z 

f 2 , . -2 

/ x da - A • x 

f xzda - A • x • z 


4A 2 j/ ^ds/t 
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1^ is the moment of inertia of A taken about an x axis passing 
through the centroid 

Iz z is the moment of inertia of A taken about a z axis passing 
through the centroid 

^2 is the product of inertia of A with respect to the centroid 
J is the tortional stiffness constant 
A^ area enclosed by the cross section 

ds incremental distance along the box contour 
t thickness of the shear resisting material associated with ds 
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2.7.3 STATIC STRENGTH ANALYSIS . The static strength analysis consists of 
calculating margins of safety for a number of failure modes. The failure modes are 
listed in Table 2-4. The analysis procedures are taken from Reference 13 . 

Simplifying assumptions are made' regarding the buckling coefficients used for calcu- 
lating local buckling stresses. Buckling coefficients for long narrow plates with 
simply supported edges or one long edge free and the others simply supported are used. 
The general instability of panel elements is checked using wide column analysis and 
assumes simple supports at ribs or frames. For multi-spar construction using thick 
plates or sandwich elements an orthotropic panel general instability analysis is made. 
The panel is checked for general buckling due to inplane loading, the lowest buckling 
load obtained from any combination of up to five half waves in each direction is used. 
The depth of analysis is generally consistent with typical predesign stress analysis. 


Table 2-4. Panel Element Failure Modes 


FAILURE MODE 

Local Buckling 

Crippling 

Inter-Rivet Buckling 

Panel General Instability 

Wide Column Buckling. . 

General Yielding 

Maximum Shear Stress 

Von Mises Equivalent Stress 

Maximum Fiber Strain in a Laminate . . . 

Diagonal Tension 

Shear Budding 


PANEL CONSTRUCTION TYPE (See Fig. 2-3] 


1 

2 

3 

4 

5 

6 

7 

8 

9 

10 

11 

12 

• 

• 

• 

• 

• 

• 

• 

• 

• 




• 

• 
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• 
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Margins of safety for the critical failure modes for each element and load condition 
are output. The margins of safety are calculated using standard interaction formula 
and are estimates of the additional load carrying capability of the element for each 
loading condition. See Section 3 of Reference 20. 

2.7.4 FATIGUE STRESS SPECTRUM. The fatigue stress spectrum is based on the 
flight profile and load spectrum discussed in Section 2.3. It is made up of a group of 
minimum and maximum stresses and the number of applications expected during the 
design life. This spectrum is used for the fatigue analysis and the flaw growth 
analysis of the element discussed in this section. 

Minimum and maximum stresses are calculated for each subsegment of the fatigue 
spectrum (see Table 2-2). These stresses are calculated from the segment constant 
stress o c and the segment alternating stress 0 a . 


c 7 . 
mm 


% ~ Ag ‘ 


0 

a 


0 

max 


<y Q + Ag • o; 


a 


The value of 0 C and 0 a are in general different for each segment and are calculated 
by forming linear combinations of the stress due to the fatigue spectrum conditions 
(see Table 2-3). 

6 





6 


a. - E 

a] 


i = 1 


a.. 0 . 
U 1 


where i denotes the fatigue spectrum condition number and j denotes the segment 
number. The constants c^ and a^ are based on the flight profile (See Table 2-1) 
and are stored within the program. 


The ground-air-ground (G-A-G) cycle shown in Figure 2-16 is not defined in Section 2.3. 
It dominates fatigue damage and flaw growth in many areas of the structure of transport 
aircraft. This stress excursion is due in part to the difference between the grotuidborn 
load distribution and the airborn distribution, and in part to cabin pressurization. 

A G-A-G spectrum is calculated automatically within the program at each analysis 
point. The G-A-G cycle is defined as the maximum stress excursion between the peak 
inflight stress (i. e. the maximum gust occurring in that flight) and the peak/valley 
groundborn stress (i.e. the maximum taxi Ag). Several high peaks of cyclic loads, 
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such as those due to gust encounters In stormy weather, tend to occur on the same 
flight. It would therefore be conservative to use all peak loads expected in the total 
aircraft life in building the G-A-G spectrum. To avoid this over conservatism a 
f requency factor is introduced which has the effect of skipping over some of the peak 
loads. A frequency factor equal to two Is considered appropriate for transport air- 
craft. Thus, every other peak is included in the G-A-G spectrum. Frequency factor 
is a user input. 



Figure 2-16. Simplified Flight Profile 


A unique stress spectrum is generated for each structural element based on the local 
stress history and is used for the fatigue analysis and flaw growth analysis also 
reported in this section. No provision is currently available for changing the load 
profile, however additional load spectra can be incorporated into the existing program 
with very little programming effort. 

2. 7. 5 FATIGUE ANALYSIS . Fatigue damage is defined as the ratio of the number of 
applied stress cycles, n, of a given stress magnitude to the number of allowable 
stress cycles, N, of the same stress magnitude. Miner's Rule (Ref. 16) is the basis 
of fatigue damage analysis performed within the program. Under this concept i 
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fatigue damage is assumed to be linearly cumulative, and fatigue failure is assumed 
to occur when the damage summation equals unity. 


Fatigue Damage 


Fatigue Failure 



n 

m 



To facilitate the analysis,S-N curves are plotted from test data for several values of 
stress ratio, R. Allowable cycles for each subsegment are read from the curves 
as shown in Figure 2-17. 

where, S 

r = ■ m - (Also see Figure 2-16) 

max 



Figure 2-17. Fatigue Damage Determination 
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A review of previous General Dynamics programs and other sources did not produce 
nearly enough component S-N data to fill the required data bank indicated by Table 2-5. 
Much of the data reviewed was generated for specific configurations and load spectra. 
Manufacturers usually test splices and other fatigue critical details but seldom develop 
S-N curves for typical structure and spectra. Even less data is published because com- 
ponent test results are frequently considered proprietary or sensitive to a particular 
project. 

Table 2-5. Availability of Fatigue Data: 


Material 

Fabrication 

Method 

Wing i 
Box 

Fuselage 

General 

Structure 

Coupon 

S-N 

Aluminum 

Riveted 

Integral 

Bonded 

Welded 

9 

9 

O 

O 

9 

O 

o 

o 

9 

9 

O 

o 

• . 


Riveted 

o 

o 

o 



Integral 

© 

o 

© 

A 

Titanium 

Bonded 

o 

vO 

o 



Welded 

o 

o 

o 



Riveted 

o 

o 

o 


Graphite/ 

Epoxy 

Integral 

o 

o 

o 

9 


Bonded 

o 

o 

o 


Boron/ 

Riveted 

o 

o 

9 


Epoxy 

Integral 

0 

o 

0 

9 


Bonded 

o 

o 

9 



O * No data 

= Complete data 


In some cases the component data was incomplete. To facilitate extrapolation, curves 
of stress vs. stress ratio at constant cycle values were plotted from the original data. 
Expanded S-N curves were then drawn based on the extrapolated data. For the many 
cases where component data was not available, reduction factors were applied to un- 
notched coupon data for the appropriate material. A complete set of data was gener- 
ated by this method. However, S-N curves plotted from this data did not show the 
trends and consistency expected. Anomalies in the component data due to inconsistent 
test parameters were still present in the expanded S-N curves. A complete and con- 
sistent set of S-N curves for all required component types could not be obtained with 
this approach. 

Subsequently, a second method, Reference 15, for plotting S-N curves from limited 
data was employed. This procedure utilizes two Hewlfett- Packard 9100B computer pro- 
grams to curve fit and plot data. More reliance is placed on un-notehed coupon data, 
and component fatigue strength factors are plotted versus life to ensure that smooth and 
consistent S-N curves are generated. All of the S-N curves shown in Appendix A were 
generated by this method. The Reference 15 procedure is presented in Appendix B for 
convenience. 
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Data from constant life cuts of these curves is stored in the program. An interpolation 
routine is used in the program to retrieve allowable cycles from the stored data. 

This rather simple approach is widely used in fatigue life predictions of transport 
aircraft. The more severe load spectra of fighter type aircraft produce more 
significant residual stresses at points of stress concentration and may warrant 
a more sophisticated analytical treatment. 

2. 7. 6 FLAW GROWTH ANALYSIS. This section presents the flaw growth analysis 
used in the VDEP II Program. The flaw growth anal ysis predicts how a crack grows 
under the influence of a fatigue load spectrum. The growth analysis is based on an 
integration technique developed under this contract. The technique currently used is 
conservative, making no provision for growth retardation effects due to spectrum 
loading. 

Crack growth predictions are usually based on the integration of empirical growth 
rate laws. The integration techniques range from simple cycle by cycle summations 
to more sophisticated techniques involving high powered numerical methods. Cyclic 
growth rate equations are usually expressed as functions of the stress intensity 
factor range, AK , the cyclic stress ratio, R, and certain empirical constants 
determined from tests data. 

— = F (AK , R, empirical constants) 
dN 

where “ is the cyclic growth rate. 
dN 

The stress intensity factor range, AK, is a measure of the change in the intensity of 
the stress field near the tip of the crack, (see Figure 2-18). 

Ak = A. (a) A (7 -/rTa = K - K 

max min 

where, 

Act range of the remotely applied cyclic stress 
a half crack length 

X(a) correction factor which accounts for geometric effects 
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The amount of growth, Aa , which occurs during one cycle of applied stress can be 
predicted using a growth rate equation, i. e. , 

A a f ( AK, R, empirical constants) 

After a load cycle is applied the new crack length Is longer by A a. at each tip. The 
AK calculated for the next load cycle is then based on the new crack size. Crack 
growth can in this way be predicted for cyclic loading. This process is cumbersome 
if the number of load cycles is very large. 



Figure 2-18. Fatigue Crack Loading 
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Crack growth analysis of transport aircraft often involves millions of cycles of 
applied loads during service life. Computer routines written to predict flaw growth 
based on the method just described were found to be too slow to use for this program. 

A technique based on averaging the cyclic growth rates during a flight and then per- 
forming the integration on a flight by flight basis was developed. A curve of average 
flight crack growth rate da/dF versus crack length is obtained for the specific crack 
geometry and load spectrum. Such a curve is obtained by summing the growth ratesidue to 
each of the load cycles in the spectrum and then dividing by the number of flights 
represented by the spectrum. By performing this process at various crack sizes a 
curve can be constructed as shown in Figure 2-19. 



The inverse of this curve, i.e. dF/da, is integrated numerically over the desired 
crack interval. A technique based on the Erdogan growth rate equation was developed 
in order to simplify and automate the procedure needed to generate the da/dF curves. 

Erdogan growth rate equation: 


“ = CK™ 
dN max 


AK P 


(2.7-1) 


where C, m, and p are empirical constants determined from constant amplitude 
crack growth tests. 


Since 


K max 


AK 

1-R 
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da _ C(AK) 

dN m 

(l-R) 


m+p 


(2.7-2) 


The average flight growth rate at any crack size a is then, 


N, 


dF M ^ 


(AK.) 


m+p 


i=l (1-Rj) 


m 


(2.7-3) 


where, 

Ak. = X (a) Act. JWsi (2.7-4) 

i i 

R. = a /a (2.7-5) 

l min. max 
i i 

Act . = ct _ CT (2.7-6) 

l max mm. 

iv . 

N g = Total number of cycles contained in the spectrum. 

M = Number of flights represented by the spectrum. 

Rewriting Equation 2. 7-3 , assuming the change in a to be small in N g cycles. 


da 

dF 



_1_ 

M 


N /A m+p 

E 


1=1 


m 


Now define CT as 



(Act.) 


m+p -i 


1 

m+p 


(l-Rj)' 


m 


(2.7-7) 


(2.7-8) 


The application: of N/M cycles of ct with an R value of 0, will produce the growth 
which is equivalent to the growth caused by the actual load cycles of an average flight. 
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da 

dF 



X(a) vJW a 


m+p 


(2.7-9) 


The number of flights required to grow a crack from a^ to a. can be found by 

solving the differential equation. Equation 2. 7-9 is separable and can be written in 
integral form as shown by Equation 2.7-10. 


AF = 


M 


x _ m+p 
N ca 
s 



a. 

i 


da 

[X (a) N /Fa] in+P 


(2.7-10) 


The closed form solution to Equation 2. 7-10 can become exceedingly complex and 
perhaps impossible depending on the form of X (a). The simplest form of X (a) which 
is useful leads to the following choice . 

X(a) = X Q +X.a (2.7-11) 

> 

Even with this simple form, the closed form of the solution is not practical. Numeri 
cal integration techniques must be applied to find the solution. Since X (a) has been 
chosen to be linear in a, the following form of Equation 2. 7-9 can be written. 


AF = 


M 


/ 


i+1 


da 


N C[ctX (a*)] m+P a. (JTa.) 


m+p 


(2.7-12) 


where 


a. 5 a* £ a. , 
i i+I 


Performing the integration in Equation 2.7-12 leads to the following solution. 


AF = 


2M 


N -C [aX (a*)] m+p * (m+p-2J 
s 


(2-m-p)/2 (2-m-p)/2 

a. —a, 

i i+I 


(2.7-13) 


where 


m+p/2 
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The error of integration introduced is dependent on the selection of a*. If a* is 
chosen to be equal to a. then the error is 


E = 



X(a*) 

m+p 

- 1 


*(a.) 



(2.7-14) 


Since (a. ^ a* ^ a. ) the maximum possible error is introduced when a* = a. 

1 i+I 7 l+l 


E 

max 



r X(a. ' 

m+p 


1+1) 

- l 


>(a.) 



(2.7-15) 


The greatest potential error can be limited to some prescribed error E* > 0 by 
proper selection of a 

i+I 



X + X • a 

m+p 


o 1 i+1 

- 1 


X + X • a. 

L O 1 1 J 



< E* 


(2.7-16) 


An operator H is introduced so that the absolute value bars can be removed. 


X + X • a. , 

1+1 < (1 + HE *)* 7 m+ ^ 
X + X • a. ' ' 

o 1 1 


(2.7-17) 


Where, 


H = 1 for X > 0 
= -1 for X 1 < 0 


Solving Equation 2. 7-17 for a results in the following 


a. 


i+1 


1/ (m+p) f 1 X o 

= (he* + d -l-r +a iJ-r 


(2.7-18) 


where X^ / 0 
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For the case where X^ = 0 the value of a^ may be set to the largest value of a 
for which the linear approximation to X is valid, in which case the integration is exact. 

The following equation is used to supply the (X q Ai) term for values of X supplied at 
two crack sizes. 



X.-a 


i±i 




(2. 7-19) 


where (a , a ) is the interval over which the X function is to be linearized and 
V j j+1 

X andX are the values of X at a and a respectively. 

J j+1 3 J+1 


The foregoing procedure is applicable to any flaw geometry for which the stress 
intensity correction factors, X, are known. The program currently contains factors 
for a wide range of stiffened panels with through cracks. Curves for L(a) and X(a) for 
the case of a crack extending equally on both sides of a riveted stiffener (illustrated 
in Figure 2-20) are stored within the program in the form of data tables. For fur- 
there information concerning the derivation of these curves the reader is referred to 
Reference 6. 



Figure 2-20. Stiffened Panel Crack Geometry 

The program currently contains 75 sets of data for L(a) and X(a) covering a wide 
range of stiffener spacing and percent stiffening, including cases for broken stiffeners. 
Figure 2-21 presents a typical set of curves. Linear interpolation is used to deter- 
mine L(a) and X(a) curves for cases which lie between data sets. These curves are 
used for all riveted- stiffener plate combinations, (e. g. , panel types 4 through9 of 
Figure 2-3). 
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Figure 2-21. Stiffened Panel Stress Intensify Correction Factors 

For the case of integral construction (e.g. , panel types 1, 2, and 3), the panel is 
treated as a flat plate without stiffeners with a thickness equal to F (i. e. , X(a) = 1.0). 

2.7.7 RESIDUAL STRENGTH. ANALYSIS. The residual strength analysis determines 
the failing strength of a damaged panel. Damage consists of skin cracks and broken 
stiffeners. The residual strength of a damaged panel is defined as the maximum stress 
level which can be applied to the panel without the crack growing unstably to failure. 
Unstable crack growth occurs when the applied stress intensify factor, K, exceeds the 
fracture toughness of the skin material,- Kq. 

Unstable crack growth is allowed to occur at stress levels below the residual strength 
of a panel as long as the crack growth eventually arrests at a larger crack size. When- 
ever stress level of the most highly loaded stiffener exceeds the ultimate tensile 
strength of the stiffener, it fails, and the applied stress intensity factors of the skin 
are recalculated to reflect the broken stiffener. 

Figure 2-22 illustrates a typical example of the residual strength analysis procedure. 
The curves shown are generated by calculating the gross panel stress which causes 
stiffener failure arid the gross panel stress which cause unstable crack growth. Equa- 
tions 2. 7-20 and 2. 7-21 are used for these calculations. 
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(2. 7-20) 


a 

cr 

(unstable crack growth) 



X(a) -/rfa 


a 


cr (stiffener failure) 



(2.7-21) 


INITIAL DAMAGE SI UtKEH, 

SI S2 S3 



Figure 2-22. Typical Example of Residual Strength Analysis 
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SECTION 3 


PART DEFINITION MODIFICATIONS 


During the course of the NASA contract NAS 1-1 1343; methods were developed to pre- 
dict a generalized detail parts listing based on the output from the structural synthesis 
analysis. The resultant parts listing provides a complete breakdown of an airframe 
structure into its most elementary components. This listing of detail parts then repre 
sents the basis of a weight and manufacturing cost analysis procedure. As part of this 
procedure each detail part is looked at individually and analyzed in terms of its weight, 
cost to manufacture, and cost of material. 

A sequence of assembling the detail parts is also modeled. The weight and costs 
corresponding to each step of the assembly process are included in the computations. 
Thus, the weight and costs of the complete airframe structure can be obtained by 
summing the values for the individual parts, and adding in the elements associated 
with the assembly. 

As part of the computer program development tasks a more comprehensive structural 
synthesis routine, APAS, was substituted for the existing lifting surface synthesis, 
BOXSIZ. This conversion required changes to the program interfaces with respect to 
both the preceeding vehicle synthesis routines and the succeeding parts definition 
routines. As part of this conversion task an effort was made to increase the degree 
of geometrical consistancy for the components being analyzed in the parts definition 
routines. 

The structural synthesis, both the BOXSIZ and the APAS versions, contain a number 
of alternative user-selected modes of construction for each of the major components 
(cover panels, spars, and ribs). Previously, only a single mode of construction 
existed for each of the components in the parts definition portion of the analysis. As 
a result of current efforts several additional structural arrangements were made 
available at the parts definition level. A summary of the various modes of construc- 
tion currently available is presented in Figure 3-1. The broken lines indicate those 
areas of analysis performed as part of this study. 
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Figure 3-1. Summary of Structural Synthesis and Parts Definition Configurations 


COMPONENT 

STRUCTURAL SYNTHESIS 

PARTS DEFINITION 


(APAS) CONFIGURATION 

CONFIGURATION 


Panels Stiffened Plate 

Metal Faced Sandwich 

Composite Faced Sandwich 

Integral Blade 

Integral Zee 

Integral Tee ' 

Joined Zee 

Joined Jay 

Extruded (Joined) Closed Hat 
Formed (Joined) Closed Hat 
Extruded (Joined) Open Hat 
Formed (Joined) Open Hat 


•-Unstiffened Plate 


►Integral Blade 
►Integral Zee 
„ Integral Tee 


Spars Integral Web 

Build Up Web 
Corrugated Web 
Integral Truss 
Built Up Truss 



"Integral Web 
. Built Up Web 


Ribs 


Integral Web 
Built Up Web 
Corrugated Web 
Integral Truss 
Built Up Truss 



Built Up Truss 


Portions of Analysis Opened Up as Part of this Study “ “ 
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3 . 1 SKIN PANE L PA RTS DE FINITION 

As part of the conversion from the BOXSIZ synthesis driver to APAS several new 
panel configurations were added to the parts definition routines. A summary of the 
panel concepts currently available in the parts definition is presented in Figure 3-2. 
Each of the lifting surface cover panels is assumed to be a complete (integral) assembly 
by itself. There are no other detail parts associated with the panels at this stage. 
Subassembly operations accounted for in the panel analysis are the spanwise and 
chordwise panel splices. Assembly of the complete box structure, as illustrated in 
Figure 3-3, is accounted for in later portions of the parts definition process. 


€ 


TYPE 1 UN STIFFENED PLATE 




TYPE 2 INTEGRAL BLADE 




TYPE 3 INTEGRAL ZEE 



T3 

- - Bl— » 

1 -I 33 !- 

J 

. ¥ J 

l t; J 


* 


TYPE 4 INTEGRAL TEE 


T3 BX — — — | — j B3 |— 

jmTiXjgL 


FIGURE 3-2. LIFTING SURFACE COVER PANEL OPTIONS 
FOR THE PARTS DEFINITION 
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SPAR CAP COVER PANEL 



FIGURE 3-3. TYPICAL MODE OF ATTACHMENT FOR 
A LIFTING SURFACE COVER PANEL. 


Input to the panel parts definition routine includes variables from both the lifting sur- 
face geometry and the structural synthesis routine. Input from the structural synthesis 
process is comprised of panel cross section dimensions of each control station. Addi- 
tional input from the geometry routine includes the number and location of the control 
stations, and the number of control station nodes. User supplied input accounts for 
specifying a panel type with an option for also defining a maximum allowable panel 
length. 

As part of the parts definition process the semi-span is divided into a number of con- 
stant length panels. In the absence of a user specified panel length, a maximum length 
of 10.06 m is assumed. Panel widths are assumed to be defined by the node locations, 
which in turn are defined in the lifting surface geometry routine. A constant number of 
panels are assumed across the surface in both a spanwise and a chordwise direction. 

Panel cross section dimensional data is computed at the actual panel endpoints by a 
linear interpolation of the geometry at adjacent control stations. The reference geo- 
metry of a typical panel in terms of the program Fortran variables is illustrated in 
Figure 3-4. Corresponding detailed dimensions for the various panel types are pre- 
sented in Figure 3-2. The panels are assumed to be spliced, using an overlap joint, 
in both the spanwise and chordwise directions, to form a complete cover ready for 
attachment to a spar and rib box structure. 
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Panel weights are computed at three different levels of consideration. A theoretical 
weight, OPWT, is computed directly from the output of the structural synthesis 
routines. The theoretical weight is the weight of the basic, idealized structural 
element. It represents an optimum value that is based on geometry of a component 
sized simply for load carrying capability. Real world manufacturing and assembly 
constraints are not considered. Typical features not accounted for are: flanges to 
serve as attachment points, clearance allowances, material widths for edge distance 
requirements, joint load path continuity, etc. 



PSTA (K+l) 


FIGURE 3-4. TYPICAL LIFTING SURFACE PANEL ARRANGEMENT 
WITH CORRESPONDING FORTRAN VARIABLES. 
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The actual weight, ACWT, reflects the actual weight of the finished part. It is com- 
puted based on the actual geometry of the finished part, and accounts for all design, 
manufacturing, and assembly considerations that would normally go into producing 
a real part. The material purchase weight, MAWT, is the weight of raw material 
stock that must be purchased in order to be able to manufacture a part of actual weight, 
ACWT. Calculation of the material purchase weight uses the same terms as the actual 
weight but includes allowances for material removed during manufacturing. Losses 
occurring between the time the material is purchased and its being utilized to produce 
a useful part (including losses due to design changes, part duplication, spoilage, waste, 
overbuy, etc.) are accounted for in the cost analysis portion of the program. Actual 
manufacturing operations which result in the loss of material include the initial mater- 
ial cut off from the raw stock, initial cutting to size, trimming, milling, turning, 
drilling, etc. Figure 3-5 illustrates the difference due to fabrication in actual and 
material purchase weights for the cover panel arrangements which are available. , 


The general form of the equation used to compute cover panel weight is: 

weight = panel length * (area + area + 4 * area ) * density 
6 12 12 

where area and, area are the cross sectional areas at the panel ends, and area 

is the cross sectional area at the panel midpoint. The panel midpoint dimensional 
data needed to compute area^ is obtained by a linear interpolation of dimensional 
data at the panel ends. 


The optimum weight for each panel is computed utilizing the dimensions output from 
the structural synthesis routines directly. The actual weight computation utilizes the 
dimensions output by the structural synthesis with extra allowances made for manu- 
facturing and assembly clearance requirements, edge distance requirements, flange 
widths for attachment, fillets, standard and minimum gages, etc. For the current 
cases involving integrally machined cover panels a minimum gage of . 081 cm is 
assumed for the outer skin portions of the panel which are to be riveted during 
assembly. Minimum gage for the remaining portions of the panel are specified by the 
user and are accounted for in the structural synthesis. For the unstiffened plate 
configuration the thickness of the skin is increased to the nearest standard material 
stock gage. Because the panels are integrally machined, there are no stiffener 
attachment flange width requirements or fastener edge distance requirements during 
assembly to consider. Fillets with a. 635 cm radius are added to the stiffener corners. 

Material purchase weights are computed based on an assumed raw material stock 
form of either flat plate or extruded flat plate. Constant dimensions are assumed 
for the raw material stock. The actual cover panel is assumed to be machined from 
the raw material stock with skin thickness and stiffener dimensions corresponding to 
the varied load conditions occurring on different parts of the surface. A constant 
(linear) rate of taper of panel dimensions is assumed between panel endpoints. Weight 
of the purchased raw material is computed by using the dimensions of the center, 
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FIGURE 3-5. PURCHASED MATERIAL FORMS FOR 
LIFTING SURFACE COVER PANEL. 
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root panel, and adding an allowance to account for machining losses. A single size of 
raw material stock is assumed for each surface and all panels for that surface are 
assumed to be machined from this. The assumed raw material forms are illustrated 
in Figure 3-5. 

The next step in the analysis procedure is a manufacturing cost computation. For 
each detail part defined by the parts definition process, a corresponding sequence of 
manufacturing operations is automatically specified. This list represents those pro- 
cesses which are required to produce the part in the shop. Each process is repre- 
sented by an equation from which a calculation is made to determine the number of 
labor hours required for that process during production of the part. These hours 
plus corresponding labor rates, overhead rates, and efficiencies form the basis of 
determining manufacturing cost. 

No additions were made to this portion of the program during the current effort. 
However, the existing analysis encompasses the fabrication of integral tee cover panel 
configurations, and is reasonably accurate for interim use with the remaining three 
panel configurations. It is recommended that future refinements to the program 
include not only the addition of new structural configuration alternatives, but also an 
updated and expanded process listing for the concepts currently available. 
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3. 2 SPAR PART DEFINITION 


As part of the conversion from BOXSIZ to the APAS synthesis driver, several new 
spar configurations were added to the parts definition routines. These configurations 
include integrally machined and builtup spar arrangements, each with either an angle 
or a tee cap. A summary of spar configurations currently available in the parts 
definition routine is presented in Figure 3-6. 

Input to the spar parts definition routines include variables from both the lifting sur- 
face geometry and the structural synthesis routine. Input from the structural synthesis 
process is comprised of spar cap and web cross section dimensions at each control 
station. Additional input from the geometry routine includes basic surface geometry, 
the number and location of the control stations, the number of ribs associated with 
the surface. User supplied input encompasses the rib and spar types, and the total 
number of spars. A minimum of two spars, a front and rear, is assumed in the 
absence of an input. A maximum of five spars are allowed. 

The program analyses each spar of a surface individually. The true length of the 
spar is computed by assuming a structural box extending from airplane centerline to 
the surface tip, minus 30. 48 cm; allowed for attachment of a tip cap. No separate 
carry through structure inside the fuselage shell is accounted for. The spars are 
assumbed to run along constant percent chord lines, and are divided into a number of 
segments if necessary to attain the full spar length. Maximum segment length may be 
input by the user. In the absence of an input integrally machined spars are assumed 
to have a maximum segment length of 4.115 m and builtup spar segments a maximum 
length of 8. 230 m. 

Each spar segment is then analyzed separately. Basic elements comprising the spar 
segment include the caps, web, web stiffeners, and rib stiffeners. Additional spar 
parts include splice plates for splicing segments, fittings for attachment of leading 
and trailing edge components, and assembly fasteners. 

The structural synthesis process provides basic dimensions for caps, web, and web 
stiffeners at each control station. Dimensions for the spar segment cap and web 
ends are derived by interpolation between adjacent control stations. A constant rate 
of taper of cap and web dimensions is assumed between segment endpoints. Web 
stiffener dimensions are assumed constant along the length of the segment, and equal 
to those dimensions arrived at by interpolation at the inboard end of the segment. 

Optimum weights for the caps, webs, and web stiffeners are computed utilizing inter- 
polated structural synthesis dimensions. The general form of equation is the same as 
that used for the cover panels. The program next adjusts these dimensions to account 
for actual manufacturing considerations. In particular, flange sizes are checked to 
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Integral Angle/Integral Blade 


Integral Tee/ Integral Blade 


Joined Angle/ Joined Angle 


Joined Tee/ Joined Angle 


Figure 3-6. Summary of the Spar Configurations Currently 
Available in the Parts Definition Routines 
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see that they are large enough to allow for proper clearances, edge distances, fasten- 
er spacing, etc. The thicknesses of components fabricated from sheet stock are re- 
adjusted to the next higher standard material gage, joint overlaps are checked, etc. 
Actual weights are then computed using the revised dimensions. 

Material weight for the integrally machined spar segment are computed by assuming 
a solid bar of material with dimensions slightly greater than the maximum spar exter- 
nal dimensions to account for cutoff and machining of the segment. The builtup spar is 
assumed to be comprised of extruded caps, sheet webs, and bentup sheet web stiffeners. 
Material weights for each of these elements are computed using the actual weight com- 
putation dimensions plus allowances for machining, cutoff, etc. 

One spar web splice plate is assumed for each spar segment except the last one. 

Splice plates are sized to fit on the spar web inside the spar cap flanges. The width 
is derived by assuming a total of four fastener rows across the splice with appropriate 
clearances, spacing, and edge distances. Thickness is set equal to the spar web thick- 
ness. The optimum weight is set equal to the actual weight in this case since the plate 
is sized initially by functional logic and not the structural synthesis. Material weight 
is computed by assuming fabrication from flat sheet stock with a thickness equal to 
the next higher standard sheet gage. Extra length and width of 5.08 cm are added 
to the actual dimensions to account for initial cutting size. 

Rib stiffeners are specified along the spar web at each rib attachment point if the local 
web stiffener spacing is greater than 30. 48 cm.. These stiffeners serve as attachment 
clips for the rib webs and are not utilized if truss type ribs are present. One set of 
rib stiffeners are specified for external spars, and two sets (one on each side of the 
spar web) for interior spars. 

Rib stiffeners are assumed to be bentup angles with a thickness and riser height equal 
to those of the web stiffeners. They are sized for an attachment flange width carrying 
a single row of fasteners. Length is set equal to the web height. Optimum weight is 
again set equal to the actual weight. Material weight is computed by assuming fabri- 
cation from the next higher standard sheet gage. The length and width are increased 
by 5. 08 cm and 2.54 cm respectively, for the material weight computation. 

Fittings are specified on the exterior spars for attachment of the leading and trailing 
edge elements. These fittings attach to spar web stiffeners directly and to the adja- 
cent skin panels through a clip and doubler arrangement. Each machined fitting has 
associated with it two clips and two doublers. 

Fittings are sized on the basis of the local spar height. A generalized fitting design 
is assumed and the actual weight is computed by deriving dimensions in proportion to 
the spar height. Appropriate dimensions and attachment flanges are checked for 
fastener allowances and readjusted if necessary. . Optimum weight is again set equal 
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to the actual weight. Material weight is computed by assuming the fitting to be machined 
from a block of thick plate. The dimensions of the block are assumed to be those of the 
maximum fitting dimensions in each direction with an additional (5.08 cm on the 
length, 2. 54 cm on the width, and a 1. 27 cm on the thickness. 

The two doublers associated with each fitting are assumed to be cut from 160 cm 
sheet. The length and width are computed as a function of the fitting size with minimum 
values set for four rows of four fasteners each. Optimum weight is set equal to actual 
weight. Material weight is computed with an additional 5. 08 cm added to the length 
and width dimensions. 

The two clips associated with each fitting are assumed to be short pieces cut from a 
tee shaped extrusion. The length and riser height are computed as a function of fitting 
size. The base flange is sized to carry a single row of fasteners along each side of 
the riser. The thickness of both the base flange and the riser is assumed to be .127 cm. 
Optimum weight is set equal to actual weight. Material weight is computed by assum- 
ing 5. 08 cm of additional length plus an allowance for the material removed from the 
riser. 

The total number of assembly fasteners required is computed by summing the number 
of fasteners necessary to assemble each detail part as it is looked at individually. 

The following assembly operations are accounted for: attachment of the segment spar 
caps to the web, splicing the webs to attain a full length spar, attachment of the web 
and rib stiffeners to the web, attachment of the doublers and clips to the skin panels, 
and attachment of the fittings. The appropriate assembly operations are bypassed 
for the case of an integrally machined spar. 

No changes were made to the manufacturi ng process listings or to the manufacturing 
cost analysis portion of the program during this effort. Spar detail parts are assumed 
to use those process listings previously specified for similar parts. It is recommended 
that future refinements include the addition of new structural configuration alternatives, 
but also an updated and expanded process listing for the concepts currently available. 



Figure 3-7. Assumed Arrangement of a Builtup Truss Type of Rib 
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3.3. RIB PART DE FINITIO N 


The existing parts definition routine for ribs was converted from the BOXSIZ synthesis 
driver to APAS directly without any additions or modifications to the mathematical 
model. A single rib configuration, a builtup truss, is available in the parts definition. 
All other rib types available in the structural synthesis revert to a builtup truss rib 
when the program reaches the parts definition process. 

The assumed arrangement of the rib is illustrated in Figure 3-7. Detail parts include 
the upper and lower (or left and right) caps, 45 degree angle braces, right angle braces, 
skin panel attachment clips, and assembly fasteners. The rib cap is assumed to be 
comprised of an extruded modified jay section with cutouts to allow for passage of the 
skin panel spanwise stiffeners. The rib caps are attached directly to the skin panels 
along the length of the caps except where spanwise skin panel splices occur. Here, 
clips are specified for attachment of the rib cap to the skin panel. The clips are 
assumed to be short pieces of extruded angle, and the number of clips associated 
with each rib is set equal to the number of spanwise panel splices. 

The number of angle braces is computed by assuming a brace angle of 45 degrees 
with a single right angle brace between adjacent 45 degree angle braces. The braces 
are assumed to be extrusions with a cruciform section. Fasteners are required for 
attachment of the braces and clips to the rib cap, and attachment of the rib cap to the 
skin panels and spars. Aluminum, steel, or titanium fasteners are available in the 
program. 

Input to the rib parts definition routine includes variables from both the lifting surface 
geometry and structural synthesis routines. The required geometry data includes the 
number and location of the control stations, the number of control station nodes, the 
number of ribs (up to a maximum of 100 are allowed per semi-span), rib spacing, 
and surface span. The structural synthesis provides at each control station the rib cap 
area and length, brace area and weight, and the structural box average height. Panel 
cross section dimensions (and panel type) are also required at each control station to 
determine the rib cap cutout size for clearance of panel stiffeners. 

Ribs are sized by the structural synthesis routine at each control station. The rib 
parts definition routine derives an actual spanwise rib location based on the specified 
rib spacing. The dimensions of the actual ribs are obtained by interpolating dimensions 
between the control stations. Each actual rib is then analyzed individually. No changes 
were made to the manufacturing process listing or to the basic rib design during this 
study. 
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SECTION 4 


CONCLUSIONS 


As a result of the work performed under this study the following conclusions are 
drawn: 

1. The overall objective of this program was to provide a user oriented computer 
program capable of assessing the impact of fatigue and fracture criteria on weight 
and cost of transport aircraft. VDEP-II meets this objective. 

2. The direct approaches adopted in the fatigue, flaw growth, and residual strengths 
analysis routines, although not specified in the work statement, were necessary to 
meet the overall objective. Allowable stress nomograph approaches were initially 
considered but they required extensive analyses external to the program for each 
selected criteria (i. e. , crack size, inspection interval, load spectrum) . Such a 
program would not be user oriented and would not provide the required accuracy, 
versatility and future growth capability. 

3. The APAS structural synthesis module developed under this study is a powerful 
tool for evaluating all types of design criteria. It enables the user to rapidly per- 
form preliminary design trade studies involving static strength factors as well as 
fatigue and fracture criteria. 

4. The sparsity of S-N data for structural components is a severe restriction on 
conventional fatigue analysis performed during preliminary design. Considering the 
wide range of parameters that effect fatigue life, it is unlikely that an adequate data 
bank for component construction types will ever be available. Unnotched coupon data 
covering the required range of cycles and stress ratios is available for most aero- 
space structural materials. Therefore, analysis methods which rely more heavily 
on this type of data should be developed. 

5. The applicability of linear elastic fracture mechanics to composite materials is 
not clearly established. Crack growth has seldom been observed. No data on growth 
rates or predictive methods was available for this program. Damage tolerance has 
been investigated for various types of construction. From the data reviewed it appears 
reasonable to treat unstiffened laminates of graphite/ epoxy by conventional fracture 
mechanics methods. This can be done with the present program by inputting the 
proper K c for each particular lay-up. 
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SECTION 5 


RECOMMENDATIONS 


During the course of the study a number of problems surfaced. These involve program 
capability, analysis methods, and data input and storage requirements. The following 
recommendations are offered to correct or alleviate these problems. 

1. The method of Appendix B, used to plot S-N Curves, should be built into the 
fatigue analysis subprogram, PRODAM. This would reduce fatigue data input and 
storage requirements, and speed up computer run time. 

2. PRODAM should be upgraded to handle biaxial stresses. At present, biaxial 
stresses generated at each section cut are converted to principal stresses before they 
are sent to PRODAM where they are treated as collinear uniaxial stresses. A more 
rational approach using octahedral stresses has been developed (Ref. 12 ) and can be 
incorporated directly in PRODAM. 

3. An input loads subprogram should be added. The present requirement to exter- 
nally generate shears, moments, and torques for six fatigue conditions plus static 
strength conditions for each structural component is very time consuming and is re- 
garded as the major hurdle for potential users of VDEP-H to overcome. Input for 
this subprogram would consist of airplane geometry, weight distributions, and aero- 
dynamic parameters. Most of this data is available in other parts of the program. 

4. The data storage approach used in PROGRO and RE SIDS should be modified to 
reduce core size requirements. At present the entire data bank of stress intensity 
factor coefficients for riveted, stiffened panels is held in core at once. This require- 
ment sets the core size for the entire program, and it can be reduced by improved 
data handling techniques. 

5. The flight profile required to develop the load spectra used in fatigue and flaw 
growth analysis should be input as a user option. At present a flight profile represen- 
tative of medium range operation of large commercial transports is installed in the 
program. A subprogram should be developed to generate load spectra from user in- 
puts consisting of weight, range, speed, altitude, and aerodynamic parameters. 

6. The program should be exercised to fully validate its accuracy and demonstrate 

its usefulness. VEDP-II as delivered has much greater analysis capability than originally 
planned and a more comprehensive checkout is required. The limited test cases which 
have been run as part of the program development are not sufficient to establish the 
needed user confidence and operating strategy. 

7. The program should be modified to include additional types of flaws commonly 
found in aircraft structure. As a minimum the following flaw types should be added: 
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a. Surface flaw 

b: Corner crack emanating from a rivet or bolt hole 

c. Single through crack emanating from a rivet or bolt hole 

d. Double through crack emanating from a rivet or bolt hole 

8. The program should be modified to automatically account for thickness effects 
on toughness. The program and data base should be modified so that as the thickness 
of structure is changed by the sizing routine, the proper toughness value is selected 
and used in the next iteration. 

9 . Flaw growth retardation due to overloads and spectrum effects should be accounted 
for. Several retardation models have been developed and one of these should be adopted 
for this program. 

10. Environmental effects on flaw growth should be added. Data for aluminum and 
titanium alloys in various environments (i. e. , R. T. air, JP fuel, LH^ should be col- 
lected and stored in the program. The input routine should allow the user to specify 
different environments for different segments of the flight. 

11. A direct method for calculating critical stress intensity factors, K c , in com- 
posite laminates should be added. Residual strength of flawed composite laminates 
can be considered in the present version of the program only to the extent that the 
user may input K c for a specific laminate. A test program is needed to define K c as 
a function of ply orientation for several composite materials. These functions should 
then be incorporated in the program so that residual strength of arbitrary angle- plied 
laminates can be calculated directly. 

12. Capability for analyzing multi- material construction types is needed. Mixing 
materials within a component has generally been avoided in past transport aircraft 
designs; however, this technique is known to be effective in dealing with crack growth 
and residual strength requirements. A procedure to optimize multi- material panels 
(e. g. , aluminum skin with titanium stiffeners) for static strength and then augment 
the section to meet fatigue and fracture criteria would be a unique and valuable capa- 
bility to add to VDEP-H. 
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APPENDIX A 


S-N CURVES 


All of the fatigue data stored in the structural synthesis routine, APAS, is shown as 
S-N curves in this Appendix. These curves were plotted with the Hewlett-Packard 
9100B computer program shown in Appendix B. The program requires input in the 
form of fatigue equation coefficients for each material and fatigue notch factors 
for each construction type. These are developed below. 

FATIGUE EQUATION COEFFICIENTS 

Unnotched coupon data is required to determine the coefficients, C, M, and n for 
Equation (B-l). All of the data needed for Aluminum and Titanium is given in Refer- 
ence 17 and is shown in Tables A-l and A-2. Values of C, m, and n shown in the 
tables were calculated by a second H-P 9100B computer program which is also shown 
in Appendix B. 

Table A-l. Fatigue Equation Coefficients for Unnotched 2024-T3 
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Table A-2. Fatigue Equation Coefficients for Unnotched Ti-6AL-4V 


= 1186 MN/m 2 (172 ksi) 
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Data for Graphite/Epoxy 0/+45/90 is given in Reference 21 for R = +. 1 only. The 
following method is used to expand this data proportional to the distribution of R values 
for unnotched Ti-6AL-4V. 







For R >R 


known' 


(-) -(b) 

\b < Titanium \ 'c 


S = S + 1 _ 

max known \b/Titanium ^composite 


For R < R^ 


nown 


S = S 
max known 


-(-L .(d) 

' d ^Titanium \ 7 i 


itanium \ /composite 


Tables A-3 and A-4 show the results of this data expansion and the calculated values 
of C, m, and n for the composite materials over the range of R = +1. 0 to R = -1. 0. 


Table A-3. Fatigue Equation Coefficients for Unnotched Graphite/Epoxy 0/^45/90 

F^ = 379 MN/m 2 (55 KSI) 
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Table A-4. Fatigue Equation Coefficients for Unnotched Boron/Epoxy 0/-45/90 

Ftu = 448 MN/m 2 (65 KSI) 
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Regarding the data for the composite materials (Boron/Epoxy and Graphite/Epoxy), it 
must be realized that an infinite variety of laminates is possible. Each layup will 
have different fatigue characteristics. A pseudo isotropic layup, 0/^45/90, was 
selected for this program because it is a common layup and some fatigue data was 
available. 

Figures A-2 thru A-5 show S-N curves for unnotched materials (2024-T3, Ti-6AL-4V 
Graphite/Epoxy, and Boron/Epoxy) plotted from the data in Tables A-l thru A-4. 
These curves are normalized to the material tensile ultimate strength, F . 



Figure A-2 . S-N Curves for Unnotched 2024-T3 
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Figure A-3. S-N Curves for Unnotched Ti-6A1-4V 



n; 


A 


XBttI 


‘ssaj^s uinuiixejAl pazrp3ui.ioN 
A -6 


Cycles to Failure, 


Figure A-4. S-N Curves for Unnotched Graphite/Epoxy 07*45/90 
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FATIGUE NOTCH FACTORS 


Fatigue behavior of structural components can be approximated by factoring unnotched 
fatigue data. The ratio of unnotched fatigue strength to notched fatigue strength at 
any number of cycles, N, is defined as the fatigue notch factor, Kf. The variation of 
Kf with cycles is given by Equation (B-5). The notch factors at the endurance 
limit, K^ e , and at one cycle, Kj s , are required input for the S-N plotter program. 
Figure B-2 is an estimate of the variation of Kf e with stress ratio, R, for several 
types of construction. Table A-5 summarizes all of the fatigue notch factors needed 
to generate S-N curves for the construction types required in this program. 

COMPONENT S-N CURVES 

The component S-N curves shown in Figures A-8 thru A-19 were generated with the 
Appendix B procedure utilizing the fatigue equation coefficients (C, m, n) in Tables A-l 
thru A- 4 and the fatigue notch factors (Kf g , K^ e ) in Table A-5. All curves are 
normalized to the net section static strength, S_. These curves provide fatigue data 

D 

for all of the component construction types currently allowed in the structural synthesis 
module, APAS. Values of constant life cuts thru these curves are stored in the FATTAB 
subroutine of APAS. New data may easily be installed by the user. The format is des- 
cribed on comment cards in FATTAB. 
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Table A-5 
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Figure A-6. S-N Curves for Riveted Aluminum Components 
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Figure A-7. S-N Curves for Integral Aluminum Components 
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Figure A-8. S-N Curves for Spot Welded Aluminum Components 
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Figure A-9. S-N Curves for Bonded Aluminum Components 
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Figure A-10. S-N Curves for Riveted Titanium Components 
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Figure A-ll. S-N Curves for Integral Titanium Components 
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Figure A-12. S-N Curves for Spot Welded Titanium Components 
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Figure A-13. S-N Curves for Bonded Titanium Components 
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Figure A-14. S-N Curves for Riveted Graphite/Epoxy, 0/i45/90, Components 
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Figure A-15. S-N Curves for Integral or Bonded Graphite/Epoxy, 0/-45/90, Components 



xbut 


‘ssajjs mnraxxBpi paztinuiJON 
A-20 


Cycles to Failure 


Figure A-16. S-N Curves for Riveted Boron/Epoxy, 0/-45/90, Components 
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Figure A-l 7. S-N Curves for Integral or Bonded Boron/Epoxy, 0/^45/90, Components 
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APPENDIX B 

AUTOMATED PLOTTING OF S-N 
DATA FOR STRUCTURAL COMPONENTS 


INTRODUCTION 

Fatigue analysis is frequently needed but seldom accomplished early in the design 
phase of aerospace structures. A major obstacle to such timely analysis is the spar- 
sity of component S-N data. There is no source of fatigue design data comparable to 
the static strength design data available in MIL-HDBK-5B. Fatigue life is affected by 
a wide range of parameters that include cyclic stress, mean stress, product form 
and orientation, temperature, environment, structural geometry, notch effects, 
metallurgical effects, and surface finish. An almost infinite variety of materials and 
fabrication methods would have to be tested to build an adequate S-N data bank. Since 
this is out of the question, the fatigue life of most aerospa.ce vehicles is verified by 
tests and analysis performed late in the development of the structure. This after-the- 
fact analysis cannot produce optimum structures. Tradeoff studies performed in the 
preliminary design phase may be compromized by lack of rational fatigue consider- 
ations. 

As a result of the circumstance outlined above the fatigue analyst must improvise 
with limited fatigue data in the early design phases. The usual approach is to select 
and plot the most appropriate data available from the literature and previous in-house 
programs. S-N curves are cross-plotted into constant life diagrams to facilitate inter- 
polation for additional stress ratio values and extrapolation to extend the cycle range. 
With some intuition and judgement S-N curves are then drawn andfaired-in to form 
preliminary design criteria for a particular component. 

A procedure for rapidly plotting S-N curves from limited data has been developed under 
a Convair IRAD study (Ref. 15). This method was used to generate the S-N curves for 
structural components shown in Appendix A. The procedure, employing two Hewlett- 
Packard 9100B Computer Programs, is presented here for convenience. The first H-P 
program is used to calculate curve- fit coefficients for unnotched coupon data. The 
second program is used to factor and plot the data as derived S-N curves for structural 
components. 

S-N EQUATION 

The characteristic sigmoidal shape of typical S-N data can be described with the 
following equation. 



B-l 


N= C 


(B-l) 



where: 


S is the maximum cyclic stress 
max J 

S is the net section static strength, F x SF 
s hi 

S g is the maximum cyclic stress at the endurance limit 
C, m, n, are data fit coefficients 
Solving for C, m, and n gives: 


C = 


WV 

(l-xAx/ 


In 


V 1< X f V 




- m- In 


n = 


1 - x 
1_ 

1 - x_ 


In I 


In 


VW ' 

M 

V X 3 / 

In 

[ y i< x i.-V | 

LW X e> J 



[ y 2 <X 2“ X e> J 


m = 


x. 


In 


1 - 


1 - x„ 


Where: 



In 


x_ 


x. 


In 


1 - x. 


1 -x. 


y = n 

s 

x = 


max 


x = — 

e S 


(B-2) 


(B-3) 


- (B-4) 
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Provided the static ultimate stress, S , and the endurance limit stress, S , are known, 

s e 

the coefficients C, m, and n can be found by substituting N and S values for three data 
points into Equations B-2, B-3, B-4. Equation 1 then describes a smooth sigmoidal curve 
which fits five data points : S , S and 3 intermediate points. Thus Equation B-l can be fitted 
to any known fatigue data sucti as the constant life diagrams presented in MIL-HDBK-5B. 

The following Hewlett-Packard 9100B program incorporates the above equations to 
solve for the coefficients C, m, and n. 


S-N DATA FIT PROGRAM OPERATING INSTRUCTIONS 

To execute the program to calculate C, m, and n perform the following steps in order. 


1 . 

PRESS "CLEAR" 


2. 

PRESS "END" 

■ - 

3. 

ENTER PROGRAM - SIDES A AND B OF MAGNETIC 

CARD 

4. 

TURN ON PRINT X 


5. 

PRESS "END" 


6. 

PRESS "CONT" 


7. 

ENTER STRESS RATIO R, "CONT" 


8. 

ENTER STATIC ULTIMATE STRESS S , "CONT" 
ENTER ENDURANCE LIMIT STRESS £ , "CONT" 
ENTER ANY THREE INTERMEDIATE &ATA POINTS 


9. 


10. 

AS FOLLOWS: 


a. ENTER CYCLES N , "CONT" AND STRESS S , "CONT" 

b. ENTER CYCLES N , "CONT" AND STRESS S , "CONT" 

m Ca 

c. ENTER CYCLES N , "CONT" AND STRESS S„, "CONT" 

o o 

11. THE PROGRAM WILL PRINT C, m, and n. 

12. THE PROGRAM RETURNS TO STEP 7. 

The programming code is listed on pages B-14, 15, and 16. 


NOTCH FACTOR EQUATION 

The ratio of unnotched fatigue strength to notched fatigue strength at any number of 
cycles, N, is defined as the fatigue notch factor, K^. Examination of notched and un- 
notched data indicates the variation of with N can be represented by the following 
inverse hyperbolic tangent equation. See Figure B-l. 


LnN = LnC, + C Tanh 
1 2 


2K -K -K 
f fs fe 

K. -K 
fe fs 


(B-5) 
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Where: 


K„ 


is the ratio of unnotched fatigue strength, to notched fatigue 
strength. 


K equals K at one cycle (static). 

tS . • • ' 7 

K equals K„ at endurance limit (assumed to be N = 10 ). 

fe • . • • ' f . • 



Figure B-l. Inverse Hyperbolic Tangent Function 


Equation B-5 is solved for as follows: 

Expressing the inverse hyperbolic function in logarithmic form gives 


Tanh 1 x = 1/2 In j 


Substituting into Equation B-5 gives 


LnN = LnC + C 
1 2 


Let a = — — - and simplify 


1/2 Ln 


V 


2K - 

K. - 

K \ 

1 + f 

fe 

fs 

K. 

- K 


fe 

fs 


2K - 

K r - 


1 - f 

fe 

fs 


K r - K 
fe fs 


J 


K ' K f 

LnN = LnC +a.ln I — — 

1 V K - K 

f fe 


(B-6) 
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And solving for yields 


V 


K + X 
fs fe 



1/a 


1 + 



(B-7) 


The component static strength factor, K^, is equal to the unnotched material ultimate 
strength divided by the component net section ultimate strength 


K _ = 
fs 


(F ) 

tu unnotched 


(F ) 

tu unnotched 


(Ffc.-SF) . (S ) 

m component s component 


where: 

F is the material ult. strength, 
tu 

SF is the net sect, factor, (i. e. , rivet factor) . 

S is the net sect, static strength 
s 


(B-8) 


The endurance limit factor, K , must either be determined from test data or esty 
mated when data is not available, A review of available data shows (at N = 10 cycles) 
to vary with the stress ratio, R, as shown in Figure B-2. At R equal plus one must 
equal . These trends may not hold for all materials, construction types, and foad 
spectra therefore data from components tested to realistic load spectra is preferred. 
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It also appears that and a can be defined as functions of R . The following 
equations are proposed. 

( R+7.07 \. 

: r~r " J 1 ”' 10 ih-9) 

V e 

a = - . 6R - 1. 4 For R > 0 (B-lOa) 

a = - 1.9R+ 1.4 For R £ 0 (B-lOb) 

The development of Equations B-5 thru B-lOb was facilitated by a Hewlett-Packard 
9100B Plotting Routine. Figure B-3 was obtained with this H-P program and is a 
plot of Equations B-6, B-9, B-lOa, and B-lOb. 


Cyoles. N 



Figure B-3. Fatigue Notch Factors for Riveted 
Alu min um Components 
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S-N PLOTTER OPERATING INSTRUCTIONS 


A plotting routine incorporating Equations B-l, B-7, B-9, B-lOa, and B-lOb has 
been programmed for the Hewlett-Packard 9100B. To operate the program execute 
the following steps: 

1. PRESS "CLEAR" 

2. PRESS "END" 

3. ENTER PROGRAM - SIDES A AND B OF MAGNETIC CARD 

4. TURN ON PRINT X 

5. TURN ON PLOTTER 

6. INSERT AND ADJUST PLOT PAPER (See Note 3). 

7. PRESS "END" 

8. PRESS "CONT" 

9. ENTER STRESS RATIO, R, "CONT" 

10. ENTER STATIC STRESS FACTOR, K, , "CONT" 

fs 

11. ENTER ENDURANCE STRESS FACTOR, K , "CONT" 

X v 

12. ENTER ENDURANCE STRESS, S e> "CONT" (See Note 2) 

13. ENTER COEFFICIENT, C, "CONT" (See Note 3) 

14. ENTER COEFFICIENT, m, "CONT" 

15. ENTER COEFFICIENT, n, "CONT" 

16. THE PROGRAM WILL PLOT AN S-N CURVE NORMALIZED TO NET 

SECTION STATIC STRENGTH, S . 

s 

17. THE PROGRAM WILL RETURN TO STEP 9. 


Notes: 

1. The program will plot both notched and unnotched S-N curves. For unnotched 

curves set K and K equal to one in Steps 10 and 11. 
fs fe 

2. The value of S entered in Step 12 must be normalized by dividing by S for the 
material from which C, m and n were calculated. 

3. The coefficients C, m, and n may be easily calculated with the H-P program 
presented previously. 

4. The program is set for semi-log paper with 7 cycles x 60 divisions. (K&E 
No. 466463). If different plot paper is used it will be necessary to change the 
scale factors. Referring to the program listing, the x and y scale factors are 
located at Steps 62 and 57 respectively. 

5. To correct anomalies in plotted curves, it may be necessary to recalculate^ 

C, m, and n using more judicious input stresses. Equation B-l is very sensi- 
tive to C, m, and n. 

The programming code is listed on Pages B-17, B-18 and B-l 9. 
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EXAMPLE PROBLEMS — DATA FIT PROGRAM 

Tables B-l and B-2 show the required input data and the calculated values of C, 
m and n for unnotched'2024-T3, and for typical riveted aluminum components. 
Also shown is a partial listing of the H-P printout. 


Table Brl. Fatigue Equation Coefficients for Unnotched 2024-T3 

F = 503 MN/m 2 (73 ksi) 
tu 


R 

NORMALIZED MAX. STRESS ,S ,F * 

max/ tu 

C 

m 

n 

Ro 

0 

4 

10 

io 5 

6 

10 


+1.0 










.8 

! 1.000 

.999 

.993 

.985 

.980 

17220 

.6674 

-103,7 

. 6 

i 


.997 

.945 

.897 

.856 

2611 

.1205 

-27.85 

.4 



.993 

.856 

.733 

.688 

13960 

.3137 

-5.103 

.2 



.984 

.753 

.589 

.551 

21140 

.3904 

-1.763 

0 



.956 

.670 

.493 

.453 

22370 

.4964 

-1.299 

-.2 



.908 

.603 

.432 

.401 

36060 

.8378 

-.3845 

-.5 


r 

.822 

.534 

.370 

.342 

73220 

1.551 

.2462 

1 

• 

o 

jl.000 

.692 

.462 

.308 

.281 


3.632 

1.206 

-2.0 










— CO 






.i ii 





* Ref. 17, page 3-77 


Sample H-P Printout: 
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Table B-2. Fatigue Equation Coefficients for Typical Riveted Aluminum Component 


S * 345 MN/m 2 (50 ksi) 
s 


R 

NORMALIZED MAX. STRESS*, S /S 

max s 

C 

m 

n 

0 

10 _ 

10 4 

■■ 


7 

10 

+1, 0 

1 . 

0 



. 





.8 

i 


.970 

.850 

.704 

.684 

501000 

1.437 

4.193 

.6 



.950 

.724 

.530 

.498 

108600 

1.050 

.6764 

.4 



.924 

.638 

.442 

.404 


.9027 

-.2568 

.2 



.882 

.552 

.350 

.318 

. 

1.157 

.2584 

0 



.818 

.482 

.296 

.266 

. 

1.413 

.2156 

-.2 



.720 

.420 

.250 

.224 

119860 

2.343 

.6162 

-.5 



.592 

.344 

. 200 

. 174 • 

130450 

3.539 

.5115 

ESI 


r 

.434 

.256 

.148 

.126 

219170 

6.532 

.6556 

EH 

1 . 

0 









* Ref. Design criteria for a transport aircraft fuselage. 


Sample H-P Printout: - . ' R 
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EXAMPLE PROBLEMS - S-N PLOTTER PROGRAM 

The S-N curves for unnotched 2024- T3 shown in Figure B-4 were plotted from the 
data shown in Table B-1. These curves are not factored, i. e. , and K were 
set equal to one. 

Figure B-5 shows S-N curves for typical riveted aluminum components plotted from 

the data in Table B-2. Again K„ and K were set equal to one. 

fs fe 

Figure B-6 is an example of component S-N curves derived by factoring unnotched 

data. The unnotched data for 2024- T3 in Table B-1 is used. The static strength 

factor, K, is found by Equation B-8. A rivet factor of 0. 75 is assumed, 
fs ! 

73 i 

K fs- 66 x 0.75 ~ 1,48 

; ; ■ . .■ * . ] 

! . " . 

The fatigue notch factors, Kf e , are from Figure B-2. The derived S-N curves 
of Figure B-6 agree within six percent with the data based S-N curves of Figure B-5. 
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Figure B-4. S-N Curves for Unnotched 2024-T3 
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Slop . 


S-N DATA FIT COEFFICIENTS - - C, m, n 
(Hewlett-Packard Programming Code) 


S-N Data Fit - - C, m, n 
Title - 
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S-N Data Fit — C, m, n 
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S-N PLOTTER FOR STRUCTURAL COMPONENTS 
(Hewlett-Packard Programming Code) 

S-N Plotter 
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